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Abstract 

Measurements of the QCD colour factors Ca and Cf and of the number of active 
quark flavours n/ in the process e + e~ — > hadrons at high energy are presented. They 
are based on fits of C(«g)+NLLA QCD calculations to distributions of the event 
shape observables 1 — T, C, Bt and Byy measured at centre-of-mass energies from 
14 GeV to 189 GeV. Hadronisation effects are approximated with power correction 
calculations which also depend on the QCD gauge structure. In this approach 
potential biases from hadronisation models are reduced. Our results for individually 
measured quantities obtained from 1 — T and C are 

n/ = 5.64 ± 1.35, C A = 2.88 ± 0.27 and C F = 1.45 ±0.27 

in good agreement with QCD based on the SU(3) symmetry group where rif = 5 
for the energies considered here, Ca = 3 and Cf = 4/3. From simultaneous fits of 
Ca and Cf with 1 — T and C we find 

C A = 2.84 ± 0.24 and C F = 1.29 ± 0.18 , 

which is also in good agreement with the QCD expectation. 



(Submitted to European Physical Journal C) 



1 Introduction 



The theory of strong interactions, Quantum Chromo Dynamics (QCD), gives a successful 
description of most aspects of hadronic final states in e + e~ annihilation and other high 
energy particle collision processes over a large range of centre-of-mass (cms) energies, see 
e.g. HIHTOf. In particular, the strong coupling a s , the only free parameter in the theory, 



can be measured to an accuracy of about 3% JTT | . 

Standard QCD contains three distinct colour charges which are carried by the quarks 
and obey the SU(3) symmetry group, see e.g. jl2|,|l3l . Requiring local gauge invariance 



of the theory generates eight massless vector fields, the gluons, each carrying a colour 
charge and an anticolour charge. Since the gluons are colour charged QCD has more 
fundamental vertices than QED, where the gauge boson, the photon, is electrically neu- 
tral. The fundamental vertices are i) the radiation of a gluon from a quark line (gluon 
bremsstrahlung) , ii) the conversion of a gluon into a quark-antiquark pair (quark pair 
production), iii) the conversion of a gluon into two gluons (triple gluon vertex) and iv) 
the conversion of a gluon into three gluons (quartic gluon vertex). The quartic gluon 
vertex is in e + e~ annihilation of order 0{a\) and thus inaccessible to the 0(a 2 s ) analysis 
presented here (see section ||). 

The relative weights of each of the vertices are given by the QCD colour factors CV, 
Tp and Ca corresponding to vertices i), ii) and iii), respectively |T2| . The colour factor 
Tp for vertex ii) contributes for each of rif active quark flavours, such that the product 
TpUf is the weight for vertex ii). 

The colour factors can be calculated from the generators of gauge transformations in 
the fundamental and adjoint representation once a particular symmetry group has been 
chosen. In standard QCD the symmetry group is SU(3) and the corresponding values of 
the colour factors are Cp = 4/3, Tp = 1/2 and Ca = 3. Perturbative QCD calculations, 
e.g. for cross sections, depend on the symmetry group only through the colour factors. 
It is thus straightforward to modify the prediction for a different symmetry group by 
recalculating the colour factors. 

From an experimental point of view one can ask the question if SU(3) is indeed the 
correct symmetry group for QCD. This question may be tested experimentally by allowing 
the colour factors to vary in comparisons of QCD predictions with data. Several such 
studies have been performed by the LEP experiments based on correlations of angles 



in 4-jet final states [14-20] . In these studies the colour factors enter at leading order 
(0(a 2 )) of the QCD calculations and thus higher order corrections have not been taken 
into account. Measurements of the ratios uj/Cf and Ca/Cp by this method have been 
in agreement with the QCD expectations. 

A complementary approach has been to employ distributions of event shape observ- 
ables measured at LEP 1 and 0(a^)+NLLA[] QCD predictions [^T],Q. Here the sensi- 



tivity to the colour factors also enters at 0{a 2 s ) but higher order corrections are partially 
taken care of by the inclusion of the NLLA terms. More recent works combined the two 



analysis techniques outlined above [23] and also made use of recently calculated higher 



order corrections to the 4-jet angular correlations [23|. In both cases results have been 



X NLLA stands for Next-to-Leading-Log Approximation, see section | 



consistent with QCD predictions and systematic uncertainties were reduced^ 

It is also possible to use the dependence of QCD predictions on the energy scale of 
the hard process to gain sensitivity to the colour factors [Il],|26|,|27|]. The dependence of 
QCD predictions on the energy scale enters through the running of the strong coupling 
a s which in turn is a function of the colour factors. 

A common property of all studies discussed so far is that they rely on Monte Carlo mod- 
els of the hadronisation process, typically JETSET JZJJ, HERWIG [g§ or ARIADNE 
Dependence on these hadronisation models in studies of the QCD gauge structure may 
be viewed as a disadvantage, because the models assume standard QCD with the SU(3) 
symmetry group to be valid. No attempt was made to study the effects of non-standard 
colour factors on the hadronisation corrections. However, such effects could be signifi- 
cant, because the parton shower stages in the Monte Carlo models are based on QCD 
calculations which in turn depend on the QCD gauge symmetry. 

Recently an analytic model of the hadronisation process has become available [|31"H 
commonly referred to as power corrections. It is based on the analysis of the power- 
behaviour of non-perturbative effects, i.e. the growth of non-perturbative effects with 
inverse powers of the scale of the hard process. In this model approximate predictions for 
the non-perturbative effects on mean values and distributions of some event shape observ- 
ables have been made and their colour structure is explicitly given [^-|37|. Comparison 
with data has been reasonably successful [p|,|5|,p7|,P^-pn|. 

In this study we use fits of £?(a^)+NLLA QCD predictions with power corrections 
to differential distributions of the event shape observables Thrust, C-parameter, Total 
and Wide Jet Broadening measured in e + e~ annihilation at cms energies from 14 GeV to 
189 GeV to investigate the gauge structure of QCD. Section |] introduces the observables 
and summarises the QCD predictions. Section | gives a description of the data sets used 
in the analysis and of the fit results. Finally in section |] a summary and conclusions are 
given. 



2 QCD Predictions 

2.1 Event Shape Observables 

We use the following event shape observables: 



Thrust T This observable is defined by the expression JTIJ [I 

\Pi-n\\ 



T = max 



where Pi is the momentum of reconstructed particle i in an event. The thrust axis 
Ht is the direction n for which the maximum occurs. We will use the form 1 — T here 
since its distribution is in this form more similar to those of the other observables. 



2 The error treatment of the results of p3] have been criticised pBj. 



C-parameter C The definition of this observable j|3],^| requires the introduction of 
the linearised energy-momentum tensor 

nal3 J^iipfPi )/\Pi\ a i o Q (n\ 

" = ^ i^i , a, p = 1,2, 3 (2) 

and its three eigenvalues A&, = 1,2,3. These define C via 

C = 3(A 1 A 2 + A 2 A 3 + A 3 A 1 ) . (3) 

Note that C may equivalently be calculated by C = (3/ J2i \Pi\) J2i<j \Pi\\Pj\ sin 2 (^j) 
where is the angle between particles i and j. 



Jet Broadening The definitions of the jet broadening observables [45| employ a plane 



through the origin perpendicular to the thrust axis Ht to divide the event into two 
hemispheres S\ and S 2 - The Total and the Wide Jet Broadening Bt and By/ are 
defined as 

£i = £p f ' |P * Wr| , B T = B 1 + B 2 , B w = max(B 1 ,B 2 ) . (4) 
2Ej |Pil 

For these observables complete 0(a 2 )+NLLA QCD predictions as well as power correction 
calculations for their differential distributions are available. 



2.2 Running of a. 



The running of the strong coupling a s is a direct consequence of the requirement that a 
complete QCD calculation, e.g. for a cross section R, should not depend on the choice 
of energy scale where the theory has been renormalised, formally expressed as the renor- 
malisation group equation (RGE), see e.g. JIB]: 



dR 







da s d 



R = 



(5) 



d/j, 2 \d/2 2 djj, 2 da 
The running of a s is found by solving the following differential equation where P(a s (n 2 )) 



is the beta function of QCD shown at two-loop accuracy p3 



2 da s 



0o 



11C A - 2n f 



P(a a (ji)) = -p a 2 s (n) - PkxKji) , 

17C 2 A - 5C A n f - 3C F n f 



(6) 



and Pi 



12tt "~ ''^ 24tt 2 
The coefficients Pq and Pi are independent of the renormalisation scheme. The solution 
to equation 



at this order is 



(7) 



with = fi/Q. The quantity is commonly referred to as the renormalisation scale 
factor since it allows to study the effects of variations of the renormalisation scale in 
perturbative QCD predictions. It is usually varied to assess the effect of missing higher 
order in the perturbation series. Equation ([7|) may be solved numerically to obtain a s {fi) 
as a function of a reference a s {Q). The QCD gauge structure enters via the coefficients 
fto and f3\ of the QCD /3-function. 

The perturbative evolution of a s diverges for small scales, see e.g. The location of 



the divergence, referred to as the Landau pole, defines the parameter A QCD which in O(o. s ) 
is given implicitly by a s (Q) = l/(/5 ln(<5 2 /AQ CD )). The value of A QCD is approximately 
200 MeV and indicates the scale where perturbative evolution of a s breaks down. 



2.3 Perturbative QCD Calculations 

The perturbative QCD prediction in 0(a 2 ) for a normalised differential event shape dis- 
tribution dRpT/dy of a generic observable y measured at cms energy ^Js = Q may be 



written as follows PTL 221: 



dRpT 1 da dA / / , 3^ \ _ dA dB\ a2 

"rf^ = = rf^^-"=^> ((-^oi n( ^) - -c-^) ac-^_ + «2C^) • (8) 

Here a to t is the cross section for the process e + e~ — > hadrons, dA/dy and dB/dy are the 
0(a s ) and the 0(a 2 ) coefficients, respectively, and a s = a s /(27r). The coefficients dB/dy 
are in fact a sum of three terms each proportional to a colour factor: 

^C F Ac A A^ . (9) 
dy dy dy J 

The coefficients dA/dy and dBi/dy, % = Cp,Cji,nf, are obtained by integrating the 
0(a 2 ) QCD matrix elements |46| in the MS renormalisation scheme using the program 
EVENT2 p7[| . A QCD prediction in 0(a 2 s ) is expected to be valid in a region of phase 
space where the radiation of a single hard gluon dominates (3- jet region). 

For the observables considered here calculations in NLLA have been performed [^, 




48|,|49|. The NLLA is valid in a region of phase space where multiple radiation of soft and 
collinear gluons from a system of two hard and back-to-back partons dominates (2-jet 
region). In the NLLA the cumulative normalised cross section Rnlla{u) is written as: 

CV 1 drT 

R N LLA(y) = / —T-dy' = C(a s )e G ^ (10) 
Jo a to t dy' 

with L = \n(l/y). The functions C(a s ) and G(a s ,L) are known in NLLA as: 

oo n+1 

C(a s ) = 1 + Cx&s + C 2 a 2 s and G{a„ L) = ^ £ G nm a n s L m ~ Lg x {La s ) + g 2 (La s ) . 

71=1 771=1 

(11) 

The coefficients C\, G n , G 12 , G22 and G23 are known analytically with complete colour 
structure while analytical expressions for the coefficients C 2 and G 2 i are not available. 



The combination of the 0(a 2 ) and the NLLA calculations is a procedure called match- 
ing. Several matching schemes have been used in measurements of a s , see e.g. IfjQfl . We 
choose to employ the so-called ln(i?)-matching scheme with the following implementa- 
tion |5T 



\nR PT (y) = L gi (La s ) + g 2 (La s ) - (G n L + G 12 L 2 )a s - (G 22 L 2 + G 23 L 3 )a 2 (12) 
+A(y)a s + (B(y) - l -A{y) 2 )a 2 s . 

The subtraction of the G nm coefficients takes account of the contributions contained in 
the NLLA as well as in the 0(a 2 ) calculations. The term A(y) is defined by A(y) = 
Jq((IA/ 'dy')dy' and equivalently B{y) = Jq(cLB / 'dy')dy' '. Our choice of the ln(i?)-matching 
scheme is motivated by two arguments, i) it is preferred theoretically j|8| as well as 
experimentally f5(| and ii) it depends only on G nm coefficients known analytically with 



full colour structure. 



2.4 Non-perturbative QCD Calculations: Power Corrections 

The model of Dokshitzer, Marchesini and Webber (DMW) lyj treats the effects of gluon 
radiation at low energy scales (O(Aqcd)) where simple perturbative evolution of a s breaks 
down due to the presence of the Landau pole. The model assumes that evolution of the 
strong coupling a s down to energies around and below the Landau pole is possible. The 
form of a s ([i) at such low energy scales is a priori unknown and a non-perturbative 
parameter is introduced as the Oth moment over a s (/j,): 

a = — / a s {k)dk . (13) 



/if Jo 

The quantity fii is referred to as the infrared matching scale where the non-perturbative 
and standard perturbative evolution of a s are merged, generally taken to be 2 GeV. 

For the observables considered here the power corrections to the differential distribu- 
tions have been calculated up to two loops |5B|,|5T|]. It turns out that a distribution dR/dy 



measured at cms energy y^i = Q can be described by the shifted perturbative prediction 
dR PT /dy: 

f = ^(V-PD V ) . (14) 
The factor P is universal [^T] and depends mainly on the non-perturbative parameter a^, 




the infrared matching scale fij and the cms energy Q: 

-a s (Q)-2(3 a 2 (Q) ^ + ^ + 1 j (15) 

with K = (67/18 — ti 2 /6)Ca — (5/9)n/ for the MS renormalisation scheme. The negative 
terms containing a s (Q) are necessary for a consistent merging of the strong coupling in 
the non-perturbative and the perturbative region. The quantity Ai (Milan factor) stems 
from two-loop effects and is given by [0 



M = 1 + (1.575C A - 0.104n / )/(47T^ ) 



(16) 



Its numerical value is 1.49 in standard QCD with n/ = 3 with an estimated theoretical 
uncertainty of 20% [fnfl . The factor D y in equation ([14]) is observable specific and is given 
by: 



Dl-T 

D c 



2 

3tt 

-In - + F y (y,a 8 (yQ)) 

^ y 



(17) 



The functions F y , y = Bt, Bw, with known colour structure, describe additional changes 
to the distributions referred to as squeeze [36. p2 . 



3 Analysis of the Data 
3.1 Event Shape Data 

We use all published data on distributions of our observables T, C, B T and B w which 
include experimental systematic uncertainties in their errors. The data available in the 
fits are listed in table [3] with references and the ranges which are considered in the fits. 
The fit ranges are determined by the following criteria: 

Data with ^ l fs > Mz° We choose to follow the approach of the OPAL collaboration [Q- 
[4|.[50|] because we use the same C(a^)+NLLA calculations with the In (R) -matching. 
The OPAL collaboration required the experimental and hadronisation corrections 
to be reasonably uniform and not strongly model dependent within the fit ranges. 
In addition, bins at the edges of fit ranges with large x 2 contributions were removed. 
The fit ranges of distributions of the other experiments are adjusted to match the 
ones of OPAL as closely as possible. 

Data with ^/s < M%o We use the studies of the recently reanalysed JADE data as a 
guideline However, we require in addition that i) the fit ranges should extend 



less far into the 2-jet region than those used in |50| and ii) the distance between the 
extreme 2-jet region and the fit ranges should increase with decreasing ^/s. These 
requirements reduce the sensitivity of the analysis to hadronisation corrections which 
are generally larger in the 2-jet region than in the 3-jet region. We also demand 
that no single bin at the edge of a fit range should have a large \ 2 contribution. 

Within these ranges experimental corrections for limited acceptance and resolution as well 
as non-perturbative effects as estimated with Monte Carlo models are well under control. 

3.2 Fit Procedure 

The fits are based on equation (|7|) for the running of a s , on equation ([T2] ) for the per- 
turbative QCD prediction and on equation (|14]) for the power corrections. We vary the 
strong coupling a s (Mzo) with the mass of the Z° boson as a reference scale, one colour 
factor, i.e. nj, Ca or Cf, and optionally the non-perturbative parameter «o in the fits. 



We also investigate fits where a s (Mzo), the two colour factors Ca and Cf and optionally 
«o are free parameters. Fixed parameters are always set to the values as expected in 
standard QCD, i.e. rif = 5, Ca = 3 and Cp = 4/3. To avoid unphysical results and to 
stabilise the fits the free parameters are bounded by 0.01 < a s (M z o) < 0.3, < rif < 20, 
< C A < 10, < C F < 10 and < a < 10. 

The fit procedure minimises a x 2 constructed from the difference in bin i between the 
data value di and the theoretical prediction ti divided by the total error af 



Possible correlations between bins of a given distribution or between different distributions 
are neglected. With the exception of the OPAL data for 1 — T, Bt and Bw at y/s = M z o 
the data are measured in bins which are wider than the typical experimental resolution in 
order to reduce bin-to-bin migrations. Many of the distributions used in this analysis were 
measured by the same experiments such that the presence of some correlation between 
distributions may be expected. There may also be correlations between the data points 
of any given distribution due to common systematic uncertainties. However, information 
about such correlations is not available in the references for the data (see table [l]). 

The errors on the fitted parameters are calculated in the fit procedure from the diagonal 
elements of the covariance matrix after the fit has converged ("fit error"). These errors 
contain the contributions from statistical fluctuations and systematic uncertainties quoted 
by the individual experiments. 

We follow two alternatives to study the influence of non-perturbative effects on our 
fits. Firstly, we fix the non-perturbative parameter a at a previously measured value and 
repeat the fits with a varied by its total errors. Secondly, we allow «o to vary in the fit. 
The fixed value of a is taken from |J, a = 0.473io!oIi- However, some modifications 
had to be applied. The measurement of || used an erroneous value of the Milan factor, 
M. = 1.794 instead of 1.49, both for rif = 3. From equation ( fL5|) we infer that it is 
sufficient to scale cto by the ratio of the two values for M. as a correction. In our study we 
vary the colour factor n / and it would therefore be inconsistent to use two different values 
of nf, i.e. rif = 5 in the perturbative part and rif = 3 in the non-perturbative part of the 
QCD prediction for massless quarks. To find our standard value of «o we scale it by the 
ratio of values for the Milan factor determined with rif = 3 and rif = 5, respectively. Our 
final value is «o = 0.543 ± 0.058 with symmetrised errors. 

3.3 Effects of bb events at low y/s 

The presence of events from the reaction e + e~ — > bb at low cms energies \fs can distort 
the event shape distributions, because the effects of weak decays of heavy B-hadrons on 
the topology of hadronic events cannot anymore be neglected. An additional potential 
problem arises from comparing QCD calculations based on massless quarks with data 
containing massive quarks at yfs close to the production threshold. 

At \/s <C M z o bb events constitute about 9% of the total event samples. Ideally one 
would correct the data experimentally by identifying bb events and removing them from 




(18) 



the sample. However, since we have only published event shape data without information 
on specific quark flavours we resort to a correction based on Monte Carlo simulations. 
We generate samples of 10 6 events at each y/s with the JETSET 7.4 program [2£| with 



the parameter set given in [53]. For each event shape observable we build the ratio of 
distributions calculated with u, d, s and c quark events to those calculated with all events. 
This ratio is multiplied with the bin contents of the data to obtain corrected distributions. 
This procedure is applied to all data at y/s < M%o. It was verified that the simulation 
provides an adequate description of the data at all values of y/s < M z o. 

Systematic effects due to uncertainties in the Monte Carlo parameters are expected 
to be small for the ratio except for those parameters which only affect the bb events 
in the samples. The most important such parameter is the value of e& in the Peterson 
fragmentation function [5"4"[ which controls the fragmentation of b quarks in the simula- 



tion. Threshold effects at low y/s which depend on the value of the b-quark mass in the 
simulation are found to be negligible for the fit results. 

3.4 Systematic Variations 

As systematic variations we consider the following changes in the analysis: 

Renormalisation scale The standard fits are carried out with the renormalisation scale 
parameter x M = 1. The dependence on the renormalisation scale is investigated by 
repeating the fits with x^ = 0.5 and x^ = 2.0 [50]. Deviations of the fit results w.r.t. 
the standard fits are taken as asymmetric uncertainties. 

Power corrections (qo fixed) The standard fits with fixed «o employ «o = 0.543, as 
explained above. The fits are repeated with ao varied by its error, i.e. «o = 0.485 
and «o = 0.601. Again, deviations of the fit results w.r.t. the standard fits are taken 
as asymmetric uncertainties. 

Power corrections (qo free) The standard fits use the Milan factor as given by equa- 
tion fllBT) . The fits are repeated with the Milan factor scaled by a factor of 0.8 and 
1.2, respectively, and deviations w.r.t. results of the standard fits are considered 
as asymmetric uncertainties. We also change the value of the infrared matching 
scale Hi from its standard value of 2 GeV to 1 GeV and 3 GeV, respectively, and 
repeat the fits. Deviations w.r.t. the standard results are counted as asymmetric 
uncertainties, except for «o where changing /// corresponds to a redefinition of the 
non-perturbative parameter chq- 

Fragmentation of b quarks The standard analysis is carried out with corrected data at 



y/s < M z o based on the JETSET tuning of |53| as explained in section |y| The value 
of the JETSET parameter e& is varied around its central value e& = 0.0038 ± 0.0010 
by adding or subtracting its error and the analysis including correction of the data 
at y/s < M z o is repeated. Deviations w.r.t. the standard results are considered as 
asymmetric uncertainties. 



Experimental uncertainties We change the composition of the data sets for the fits in 
two ways, i) the data measured at y/s = M z o are removed and ii) only data measured 



at y/s > M z o are used. The larger of the two deviations w.r.t. the standard result 
observed with the data set i) or ii) is considered as the systematic uncertainty from 
variations in the input data. 

The systematic uncertainties from the renormalisation scale, the power corrections with 
ao either fixed or free, the fragmentation of b quarks and from the input data are added 
in quadrature with the error from the fit to arrive at the total error. 



3.5 Fit Results 

The tables §, § and [| present the results of two-parameter fits to a s (M z o) and n/, Ca 
or Cp, respectively, with a fixed. The results of three-parameter fits to a s (M z o), Ca 
and Cp are given in table |5j. The results of three-parameter fits to a s (M z o), a and 
one of the colour factors rif, Ca or Cp are shown in tables [7| and respectively, 
while table |9] displays the results of four-parameter fits to a s (M z o), ao, Ca and Cp. The 
rows labelled "hadr." contain the quadratic sum of the uncertainties from the power 
corrections and the fragmentation of b-quarks. The power correction uncertainties are 
the dominating contribution in all cases. The tables show symmetrised errors except for 
the fit to a s (M z o) and nj where the variation of ao leads to a large difference between 
the positive and negative error. Figures |l|, |] and § give a graphic display of the results. 

We observe stable fits in all cases for the observables 1 — T and C. The fit results 
from these observables for the colour factors are generally in good agreement with the 
expectations from standard QCD within their total uncertainties. The results for a s (M z o) 
and «o are also generally consistent with previous measurements within their total uncer- 
tainties |6, nj|. The fit of n/ and a s (M z o) with 1 — T is consistent with an earlier analysis 



using similar data |2H]. The fit of a s (M z o), a , Ca and Cp with C converges but has such 
large fit errors that sensitivity to the colour factors is essentially lost; consequently we do 
not show the results. 

With the other two observables, Bp and Bw, we also find generally consistent results 
with standard QCD and previous measurements. However, fits of Bp to a s (M z o) and 
rif and fits of of Bw to a s (M z o), Ca and Cp as well as to a s (M z o) and single colour 
factors with a free are unstable and we do not show these results. Our interpretation is 
that the QCD calculations with power corrections for 1 — T and C provide an adequate 
description of the data while there are still effects in the data for Bp and Bw which are not 
well described by the C(a^)+NLLA QCD calculations with power corrections |5|,[^, ^0| , |50[| . 

The values of x 2 /d.o.f. are smaller than one in all cases, except the fit of a s (M z o) and 
Ca with Bp where x 2 /d.o.f. = 1.01. The small values of x 2 /d.o.f. are consistent with the 
good description of the data by the fitted predictions and probably indicate the presence 
of correlations between the data points used in a given fit. 

We find that in all fits with ao free the variation of the Milan factor leads to sig- 
nificantly smaller systematic uncertainties than the variation of ao in fits with ao fixed. 
The fitted values of ao vary over ranges compatible with the total error of the fixed value 
of ao when the Milan factor is changed. This is consistent with |J where the error on 
a was dominated by the variation of the Milan factor. We conclude that the ratio of 
non-perturbative and perturbative contributions is well constrained by the data such that 



changes in the Milan factor are compensated by similar changes with the opposite sign 
in ckq. As a consequence we find that the hadronisation systematic uncertainties of fits 
to the colour factors are reduced when ao is a free parameter while the fit errors increase 
due to the presence of an extra free parameter in the fit. 

We observe that the experimental uncertainties estimated by repeating the fits with 
reduced data sets are larger in the fits with «o free. This effect is likely to be a reflection 
of the increased fit errors. We choose to keep the larger uncertainties since we cannot rule 
out the possibility that there are systematic effects in the data to which the fits with «o 
free are more sensitive than the fits with «o fixed. 

In simultaneous fits with several free parameters correlations between them may influ- 
ence the results. In table ITOl we show the correlations observed in the standard fits with 



the observable 1 — T. In addition in the first row of table [10] the correlation observed 
in a fit with only a s (M z o) and a as free parameters is displayed. The presence of large 
correlations in some fits causes increased fit errors and may also contribute to the experi- 
mental uncertainties. We find a similar pattern of correlations with the other observables 
C, Bt and Bw- 



3.6 Combination of Results 

We construct combined results only from fits to 1 — T and C, because we found these to 
be stable in all cases. We choose to base the combination for individual colour factors on 
the fits with a free since these results should be less biased by input values measured 
assuming standard QCD. They also have smaller hadronisation uncertainties as discussed 
above. We build unweighted averages of the fit results for the standard analysis and for all 
systematic variations. As the fit error we choose the smaller of the two individual fit errors 
and we then construct the total error in the same way as for an individual observable. 
In this way correlations between systematic variations of the observables are taken into 
account. We find as our final results for the individually measured colour factors: 

n f = 5.64±0.79(fit) ±0.32(scale) ±0.48(had.) ±0.93(exp.) (19) 

= 5.64 ±1.35, 

C A = 2.88±0.16(fit) ±0.06(scale) ±0.10(had.) ±0.18(exp.) 

= 2.88 ±0.27 and 

C F = 1.45 ±0.21 (fit) ±0.04(scale) ±0.04(had.) ±0.16(exp.) 

= 1.45 ±0.27. 



The total uncertainties for all three colour factors are dominated by the fit error and the 
experimental uncertainties. The results are in good agreement with the expectation from 
standard QCD with five active quark flavours. The results for a s (Mzo) and cto from the 
combination are consistent with the individual fits and with previous measurements PJTIfl, 



see table 11 



In the case of rif the combined result has a slightly larger total uncertainty than the 
individual result based on the observable 1 — T shown in table ^. This is mainly due to 



the increased experimental uncertainty of the combined result which in turn is dominated 
by the experimental uncertainty of the fit with C. For the other colour factors Ca and 
Cp the total uncertainties of the combined results are somewhat smaller than the total 
uncertainties of the individual results. These observations are consistent, since we do not 
expect the total uncertainties to be greatly reduced by a combination of individual results 
due to the large correlations between the event shape distributions. 

A combined result for the results of the simultaneous fits of Ca and Cp may only be 
constructed from the fits with «o fixed, because in the fits with «o free we obtained usable 
fits only for the observable 1 — T . Applying the same procedure as outlined above yields 
the following results for combining the results of the fits with the observables 1 — T and 
C shown in table |: 

C A = 2.84±0.13(fit) ±0.06(scale) ±0.11(had.) ±0.15(exp.) (20) 

= 2.84 ±0.24 and 

C F = 1.29 ± 0.07(fit) ±0.17(had.) ± 0.02(exp.) 

= 1.29 ±0.18. 



The scale uncertainty of the combined result for Cp is smaller than 0.01 consistent with 
the small scale uncertainties of the individual results, see table |5]. The results for C A and 
Cp are in good agreement with the combined results discussed above, see equation (|19|) . 
The combined result for a s (Mzo) is in agreement with the individual results and also 



with 111, see table 11 



The averaged correlation coefficient from the fits is pc A -c F — —0.89. However, we 
observe that variations of the results for Ca and Cp are positively correlated when a is 
changed. In order to find a conservative estimate of the correlation coefficient we construct 
a covariance matrix by summing i) the covariance matrix from the fit and ii) one covariance 
matrix for each systematic uncertainty. The covariance matrices ii) are constructed from 
the systematic uncertainties, symmetrised if necessary, and pc A -c F — —0.89, except for 
the systematic uncertainty from the variation of ao where pc A -c F = 1-0 is used. The 
resulting total correlation coefficient is pc A -c F = 0.19. Figure |] presents the combined 
results which are in good agreement with standard QCD within the uncertainties. Some 
other possible gauge groups are indicated including U(l) 3 , an abelian QCD with three 
quark colours and colour neutral gluons. Our results exclude all shown alternatives to 
SU(3) as the gauge group of QCD at more than 95% confidence level. 



4 Summary and Conclusions 

We have studied fits of C(a^)±NLLA QCD predictions for distributions of the event 
shape observables 1 — T, C, Bp and By/ with power correction calculations to model 
hadronisation effects to data measured at cms energies ranging from 14 GeV to 189 GeV. 
In these fits we varied simultaneously the strong coupling a s (M z o), one of the QCD colour 
factors Ca or Cp, or the the number active quarks rif, and in some cases also a , the free 
parameter of the power correction calculations. We investigated in addition fits, where 
a s (M z o), Ca and Cp and optionally a were varied. 



We found stable fits in all cases with the observables 1 — T and C while some fits 
with Bt and in particular Byy appear unreliable. We take this as an indication that the 
current (9(a^)+NLLA QCD calculations with power corrections for Bj- and By/ describe 
the data not as well as the same calculations for 1 — T and C. 

We observed that the variation of the Milan factor M. and the infrared matching scale 
fii in the fits with «o free gave rise to smaller systematic uncertainties than the systematic 
variation of «o in the fits with «o fixed. Our conclusion is that the relative contributions 
of the perturbative (9(a^)+NLLA QCD calculations and the power correction terms to 
the total predictions are quite well constrained by the data. 

A combination of the individual results of the fits with a s (M z o), an individual colour 
factor and cto as free parameters with the observables 1 — T and C yields our final results: 

n f = 5.64 ±1.35, 

C A = 2.88 ±0.27 and 

C F = 1.45 ±0.27. 



The combination of results from simultaneous fits of Ca, Cp and a s (M z o) for ao fixed 
with the observables 1 — T and C gives 

C A = 2.84 ±0.24 and C F = 1.29 ±0.18 (21) 

with total correlation coefficient pc A -c F = 0.19. These results are in good agreement 
with the expectation from standard QCD based on the SU(3) symmetry group for e + e~ 
annihilation data at high energies, i.e. nf — 5 with Tp = 1/2, Ca = 3 and Cp = 4/3. 
There is also good agreement between the individual and the simultaneous measurements 
of Ca and Cp. All combined results for a s (M z o) and cto are in agreement with previous 
measurements. We found that the total uncertainties of the combined results are ap- 
proximately of the same size as the total uncertainties of the individual results. This is 
consistent, because the event shape distributions are correlated with each other and the 
total uncertainties are dominated by the fit errors and experimental uncertainties. 

We present our final conclusions from two points of view. Firstly, we assume that 
the power correction calculations are a good model of hadronisation effects in event shape 
distributions. In this case we have performed a complementary test of the gauge structure 
of QCD with competitive uncertainties on the measurements of the QCD colour factors 
compared to other analyses [[19| , |20| , |23| , p4| . Secondly, under the assumption that QCD with 



SU(3) group symmetry is the correct theory of strong interactions, our analysis provides 
a successful consistency check of the power correction calculations. 
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Tables 





Experiment 


1 — T 


C 


Bt 


Bw 


189 


L3 
OPAL [§ 


0.05-0.30 
0.04-0.30 


0.15-0.60 
0.18-0.60 


0.08-0.26 
0.075-0.25 


0.045-0.195 
0.05-0.20 


183 


DELPHI H 
L3|?J 
OPAL [| 


0.04-0.28 
05-0 30 
0.04-0.30 


0.16-0.64 
15-0 60 
0.18-0.60 


0.07-0.24 
08-0 26 
0.075-0.25 


0.05-0.20 
045-0 195 
0.05-0.20 


172 


DELPHI 
L3 H 
OPAL @ 


0.04-0.24 
0.05-0.30 
0.04-0.30 


0.16-0.64 
0.15-0.60 
0.18-0.60 


0.08-0.27 
0.08-0.26 
0.075-0.25 


0.04-0.17 
0.045-0.195 
0.05-0.20 




DELPHI [38 
L3 |7J 
OPAL [1 




04-0 24 
05-0 30 
0.04-0.30 


16-0 64 

\J . ±.\J \J . V.' t: 

15-0 60 
0.18-0.60 


08-0 27 
08-0 26 

\J • \J V7 « ^ V7 

0.075-0.25 


04-0 1 7 
045-0 195 
0.05-0.20 


133 


ALEPH [551 
DELPHI [38 
L3 |3 
OPAL [| 




04-0 30 
0.04-0.24 
05-0 25 
0.04-0.30 


0.16-0.64 
15-0 64 
0.18-0.60 


0.08-0.27 
08-0 26 
0.075-0.25 


0.04-0.17 
045-0 195 

V 7 • V 7 J- '..J V 7 • _L '7 '.7 

0.05-0.20 


91 


ALEf 
DELI 

L3 m 

OPAI 
SLD 


>H 

m 

i 

60 


15 
11 

30. 


37 

51 




06-0 30 
0.06-0.30 
0.065-0.29 
0.06-0.30 
0.06-0.26 


20-0 64 

V 7 • V 7 V 7 • V / -L 

0.20-0.64 
0.22-0.64 
0.20-0.64 
0.24-0.64 


0.09-0.24 

0.09-0.23 
0.08-0.26 


0.07-0.17 

0.07-0.17 
0.08-0.20 


55 


AMY [0 


0.10-0.30 








44 


JADE @,|6 

TASSO [BS 


1 


0.08-0.30 
0.08-0.28 


0.24-0.58 


0.10-0.24 


0.08-0.18 


35 


JADE |,1 
TASSO §2| 


0.08-0.30 
0.08-0.28 


0.24-0.58 


0.10-0.24 


0.08-0.18 


29 


hrs [g 

MARKII [|| 


0.10-0.30 
0.10-0.30 








22 


TASSO |! 


!l 


0.10-0.28 








14 


TASSO |2| 


0.12-0.28 









Table 1: The sources of the data and the fit ranges for the observables 1 — T, C, B T 
and B\y are shown. The cms energy a/s at which the experiments analysed their data is 
given in GeV. 





1 — T 


c 


Bw 




n f 


a s (M z o) 


n f 


a s (M z o) 


n f 


a s (M z ») 


fit result 


7.32 


0.124 


5.82 


0.114 


2.06 


0.093 


fit error 


±0.49 


±0.002 


±0.52 


±0.002 


±1.04 


±0.003 


xVd.o.f. 


130.8/236 


104.0/155 


109.0/121 


ren. scale 


±0.11 


±0.003 


±0.55 


±0.002 


±1.33 


±0.002 


hadr. 


+0.56 
-2.39 


±0.007 


±2.30 


±0.006 


±1.47 


±0.003 


exp. 


±0.42 


±0.003 


±1.95 


±0.006 


±2.89 


±0.010 


tot. error 


+0.86 
-2.48 


±0.008 


±3.14 


±0.009 


±3.66 


±0.011 



Table 2: Results are shown for fits with the observables 1 — T, C and B w to a s (M z o) 
and rif. 





1 — T 


c 


Bt 


B w 




C A 


a s {M z o) 


C A 


a s {M z o) 


C a 


a s {M z o) 


C A 


a s (M z o) 


fit result 


2.62 


0.123 


2.88 


0.113 


3.82 


0.101 


3.53 


0.094 


fit error 


±0.08 


±0.002 


±0.07 


±0.002 


±0.10 


±0.002 


±0.16 


±0.003 


xVd.o.f. 


128.8/236 


103.7/155 


139.6/138 


105.9/121 


ren. scale 


±0.03 


±0.004 


±0.08 


±0.003 


±0.21 


±0.003 


±0.21 


±0.001 


hadr. 


±0.42 


±0.006 


±0.49 


±0.008 


±0.36 


±0.004 


±0.23 


±0.003 


exp. 


±0.07 


±0.003 


±0.31 


±0.005 


±0.23 


±0.003 


±0.67 


±0.009 


tot. error 


±0.43 


±0.008 


±0.59 


±0.010 


±0.49 


±0.006 


±0.76 


±0.010 



Table 3: Results are shown for fits with the observables 1 — T, C, B T and B w to a s (M z o) 
and Ca- 





1 — T 


C 


Bt 


Bw 




c F 


a s (M z o) 


c F 


a s (M z o) 


c F 


a s (M z o) 


c F 


a s (M z o) 


fit result 


1.11 


0.128 


1.28 


0.114 


1.65 


0.100 


1.98 


0.078 


fit error 


±0.06 


±0.004 


±0.03 


±0.002 


±0.03 


±0.002 


±0.09 


±0.003 


xVd-o.f. 


133.4/236 


103.6/155 


131.3/138 


87.8/121 


ren. scale 


±0.05 


±0.004 


±0.05 


±0.003 


±0.01 


±0.003 


±0.07 


±0.001 


hadr. 


±0.24 


±0.013 


±0.22 


±0.009 


±0.21 


±0.006 


±0.24 


±0.006 


exp. 


±0.10 


±0.008 


±0.04 


±0.001 


±0.17 


±0.008 


±0.27 


±0.008 


tot. error 


±0.28 


±0.016 


±0.23 


±0.010 


±0.27 


±0.011 


±0.38 


±0.011 



Table 4: Results are shown for fits with the observables 1 — T, C, B T and B w to a s (M z o) 
and C F . 





1 — T 


c 




C A 


c F 


a s (M Z o) 


C A 


c F 


a s (M z o) 


fit result 


2.72 


1.28 


0.124 


2.96 


1.30 


0.114 


fit error 


±0.13 


±0.07 


±0.002 


±0.18 


±0.08 


±0.002 


X 2 /d.o.f. 


127.9/235 


103.6/154 


ren. scale 


±0.04 


±0.00 


±0.004 


±0.08 


±0.00 


±0.003 


hadr. 


±0.10 


±0.17 


±0.010 


±0.12 


±0.16 


±0.009 


exp. 


±0.11 


±0.13 


±0.005 


±0.40 


±0.10 


±0.004 


tot. error 


±0.20 


±0.22 


±0.012 


±0.46 


±0.21 


±0.011 










c A 


c F 


a s (M z o) 


fit result 


3.39 


1.58 


0.097 


fit error 


±0.21 


±0.06 


±0.002 


X 2 /d-o.f. 


128.0/137 


ren. scale 


±0.10 


±0.01 


±0.002 


hadr. 


±0.16 


±0.19 


±0.008 


exp. 


±0.48 


±0.02 


±0.005 


tot. error 


±0.56 


±0.20 


±0.010 



Table 5: Results are shown for fits with the observables 1 — T, C and Br to a s (M z o), 
Ca and Cp. 





1 — T 


C 




n f 


a s (M z o) 


a 


n f 


a s (M z o) 


a 


fit result 


6.39 


0.121 


0.521 


4.88 


0.111 


0.526 


fit error 


±0.79 


±0.003 


±0.015 


±1.34 


±0.004 


±0.023 


xVd.o.f. 


127.9/235 


103.4/154 




ren. scale 


±0.21 


±0.004 


±0.001 


±0.42 


±0.003 


±0.005 


hadr. 


±0.47 


±0.001 


±0.057 


±0.48 


±0.001 


±0.065 


exp. 


±0.69 


±0.004 


±0.010 


±2.55 


±0.008 


±0.038 


tot. error 


±1.17 


±0.006 


±0.060 


±2.96 


±0.010 


±0.079 




Bx 






n f 


a s (M z o) 


a 


fit result 


3.77 


0.108 


0.635 


fit error 


±1.12 


±0.003 


±0.019 


xVd.o.f. 


132.1/137 


ren. scale 


±0.46 


±0.004 


±0.006 


hadr. 


±0.33 


±0.001 


±0.089 


exp. 


±2.97 


±0.009 


±0.002 


tot. error 


±3.22 


±0.010 


±0.091 



Table 6: Results are shown for fits with the observables 1 — T, C and B T to a s (M z o), 
a and rif. 





1 — T 


C 




c A 


a s (M z o) 


a 


C A 


a s (M z o) 


a 


fit result 


2.73 


0.121 


0.528 


3.03 


0.111 


0.525 


fit error 


±0.16 


±0.003 


±0.020 


±0.25 


±0.004 


±0.029 


Y 2 /d.o.f. 


128.1/235 


103.4/154 




ren. scale 


±0.05 


±0.004 


±0.001 


±0.08 


±0.003 


±0.002 


hadr. 


±0.10 


±0.001 


±0.057 


±0.09 


±0.001 


±0.064 


exp. 


±0.14 


±0.004 


±0.012 


±0.50 


±0.008 


±0.053 


tot. error 


±0.24 


±0.006 


±0.061 


±0.57 


±0.009 


±0.089 




B F 






c A 


a s {M z o) 


a 


fit result 


3.24 


0.108 


0.626 


fit error 


±0.21 


±0.003 


±0.023 


xVd.o.f. 


131.8/137 


ren. scale 


±0.07 


±0.004 


±0.006 


hadr. 


±0.06 


±0.001 


±0.088 


exp. 


±0.60 


±0.009 


±0.086 


tot. error 


±0.64 


±0.010 


±0.126 



Table 7: Results are shown for fits with the observables 1 — T, C and B T to a s (M z o), 
ccq and Ca- 





1 — T 


C 




c F 


a s (M z o) 


a 


c F 


a s (M z o) 


a 


fit result 


1.42 


0.113 


0.478 


1.49 


0.105 


0.490 


fit error 


±0.21 


±0.009 


±0.041 


±0.31 


±0.011 


±0.074 


xVd.o.f. 


130.4/235 


103.2/154 




ren. scale 


±0.03 


±0.006 


±0.013 


±0.10 


±0.006 


±0.020 


hadr. 


±0.03 


±0.001 


±0.053 


±0.06 


±0.002 


±0.059 


exp. 


±0.21 


±0.011 


±0.022 


±0.10 


±0.004 


±0.031 


tot. error 


±0.30 


±0.015 


±0.072 


±0.35 


±0.014 


±0.103 




B F 


Bw 




c F 


a s {M z o) 


a 


c F 


a s (M z o) 


a 


fit result 


2.28 


0.085 


0.406 


3.50 


0.052 


0.311 


fit error 


±0.47 


±0.009 


±0.080 


±0.73 


±0.009 


±0.068 


xVd.o.f. 


129.4/137 


79.1/120 


ren. scale 


±0.88 


±0.029 


±0.222 


±0.44 


±0.005 


±0.031 


hadr. 


±0.32 


±0.006 


±0.056 


±0.38 


±0.004 


±0.051 


exp. 


±0.80 


±0.018 


±0.248 


±1.44 


±0.023 


±0.284 


tot. error 


±1.32 


±0.036 


±0.348 


±1.72 


±0.025 


±0.298 



Table 8: Results are shown for fits with the observables 1 — T, C, B T and B w to a s (M z o), 
a and C F . 





1 — T 




C A 


c F 


a s (M z o) 


a 


fit result 


2.68 


1.21 


0.128 


0.567 


fit error 


±0.19 


±0.26 


±0.016 


±0.093 


xVd.o.f. 


127.8/234 


ren. scale 


±0.03 


±0.02 


±0.005 


±0.006 


hadr. 


±0.16 


±0.14 


±0.010 


±0.042 


exp. 


±0.17 


±0.35 


±0.025 


±0.164 


tot. error 


±0.30 


±0.46 


±0.032 


±0.194 



Table 9: Results are shown for fits with the observables 1 — T to a s (M z o), a , Ca and 
C F . 



fit type 




a 


n f 


C A 


C F 


a s (M z o)-a 


a s 


-0.87 








a s (M z o)-n f 


a s 




0.985 






a s (M z0 )-C A 


a s 






-0.981 




a s (M z o)-C F 


a s 








-0.997 


a s (M z o)-C A -C F 


a s 
C A 






0.20 


-0.68 
-0.85 


a s (M z o)-a -n f 


a s 
a 


0.69 


0.97 
0.82 






a s (M z o)-a -C A 


Oi s 

a 


0.78 




-0.97 
-0.89 




a s (M z o)-a -C F 


Ois 

a 


0.96 






-0.998 
-0.975 


a s (M z o)-a -C A -C F 


C A 


0.988 




-0.72 
-0.75 


-0.985 
-0.970 
0.59 



Table 10: The correlations between the fit parameters for all types of standard fits to 
the 1 — T distributions are shown. 





n f 


a s (M z o) 


a 


C A 


a s (M z o) 


a 


result 


5.64 


0.116 


0.524 


2.88 


0.116 


0.526 


error 


±1.35 


±0.005 


±0.064 


±0.27 


±0.005 


±0.067 




c F 


a s (M z o) 


«o 


C A 


c F 


a s (M z o) 


result 


1.45 


0.109 


0.484 


2.84 


1.29 


0.119 


error 


±0.27 


±0.013 


±0.075 


±0.24 


±0.18 


±0.10 



Table 11: The combined results based on 1 — T and C are shown with total errors. The 
results for Ca, C f and a s (M z o) are from fits with a fixed. 
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Figure 1: The figures present results of fits to a s (Mzo) and one of the colour factors n/, 
Ca or Cf with observables as given on the vertical axis. The vertical dotted lines indicate 
the expectation from standard QCD for the colour factor. 
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Figure 2: The figure presents results from fits to a s (M z o) and the colour factors Ca and 
Cf with observables as given on the vertical axis. The error bars show total uncertainties. 
The vertical dotted lines indicate the expectations from standard QCD for the colour 
factors. 
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Figure 3: The figures present results from fits to a s (M z o), «o and one of the colour 
factors n/, Ca or Cf with observables as given on the vertical axis. The error bars show 
total uncertainties. The vertical dotted lines indicate the expectation from standard QCD 
for the colour factor. 
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Figure 4: The figure presents the combined results for the colour factors Ca and Cp 
from fits to a s (M z o), Ca and Cp based on the observables 1 — T and C. The square and 
triangle symbols indicate the expectations for Ca and Cp for different symmetry groups. 



